Semana 6 (Semestre 2)
Para esta semana, se trato de arreglar el problema de la fuerza armonica angular. Esto se consiguio modificando la direccion de la fuerza con productos cruz.
Se puede observar como las particulas de la cadena de abajo se 'doblan', a pesar que tienen total libertad de girar en el sistema. Esto se puede deber a que el potencial usado en la interaccion entre cadenas es el Lennard-Jones completo. Para esto, se intentara implementar la teoria de perturbacion de Weeks-Chandler-Andersen perturbation theory.
La cual esta basada en el potencial lennard jones, pero que esta dividida entre repulsion e atraccion. Para este experimento, trataremos usar meramente de arepulsion entre cadenas.
Debido a que esta fuerza solamente actua cuando la distancia es demasiada corta, en la simulacion se aprecia que ya no existe alguna fuerza entre particulas, lo cual.
Interacción magnética
Se implementó la fuerza magnética entre cadenas, la cual se hizo de esta manera
function dipole_magnetic_force(r, m1, m2)
μ = 4 * π * 10^(-7)
r_norm = sqrt(sum(r .^ 2))
force_direction = r / r_norm
m1_dot_r = dot(m1, r)
m2_dot_r = dot(m2, r)
m1_dot_m2 = dot(m1, m2)
prefactor = ( 3 * μ ) / ( 4 * π * r_norm ^ 3)
magnitude = prefactor * ((m1_dot_r * m2) + (m2_dot_r * m1) + (m1_dot_m2 * r) - (5 * m1_dot_r * m2_dot_r)/(r_norm .^ 2) * r)
return magnitude
end
Lo cual en un principio ya que al parecer las particulas estaban muy alejadas unas de las otras, no se percibia una interaccion magnetica, es por eso que al cambiar las distancias se obtuvo esto.
Donde se alcanza a observar que una de las particulas se junta con otra, esto puede ser esperado debido a que los momentos magneticos estan apuntando en la misma direccion.
Pero, que pasara si los momentos estan apuntando en la direccion contraria?.
